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Abstract 



v^ ' The 'renormalon' or 'dispersive' method for estimating non-perturbative 



corrections to QCD observables is reviewed. The corrections are power- 
suppressed, i.e. of the form A/Q^ where Q is the hard process momentum 
l^^_^ scale. The renormalon method exploits the connection between diver- 

^j< . gences of the QCD perturbation series and low-momentum dynamics to 

predict the power, p. The further assumption of an approximately uni- 
versal low-energy effective strong coupling leads to relationships between 
the coefficients A for different observables. Results on 1/Q^ corrections 
O ' to deep inelastic structure functions and 1/Q corrections to event shapes 

D , are presented and compared with experiment. Shape variables that could 

tS^ • be free of 1/Q and as{Q^)/Q corrections are suggested. 



1 Introduction 

The question of how best to compare QCD predictions with strong interaction 
data is an old one. The problem is that we are usually comparing predictions 
about quarks and gluons - the so-called parton level - with data on mesons 
and baryons - the hadron level. We must then choose carefully those quantities 
for which the differences between partons and hadrons are not expected to be 
important. 

Our criteria for choosing what to compare are in fact all based on parton- 
level perturbation theory, so the argument looks in danger of being circular. 
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However, we believe that the conversion of partons into hadrons - hadroniza- 
tion - is a long-distance process. So we can select those observables that are 
insensitive to long-distance physics (infrared-safe quantities) or in which long- 
distance contributions can be separated out (factorizable quantities). These 
properties can be established entirely within perturbation theory, and it seems 
reasonable to expect that they are also true in the real world. 

Increasingly, though, we need to know more quantitatively which are the 
best QCD observables, with the least sensitivity to long-distance physics. One 
way to study this is to make a model for the non-perturbative hadronization 
process and to look at the size of the resulting hadronization corrections. This 
has been a very fruitful approach, and is still the standard procedure of most 
experiments, but it does have its limitations. [] It doesn't apply to long-distance 
physics other than hadronization, for example the parton distributions inside 
colliding hadrons. It's a 'black-box' procedure that often doesn't provide much 
insight, and has many adjustable parameters that may conceal wrong physics 
assumptions. Perhaps most seriously of all, it doesn't really apply in the circum- 
stances in which it is normally used. By that I mean that the 'parton level' to 
which the hadronization model is applied is a somewhat mutilated version of the 
final states that appear in a perturbative calculation. It contains approximate 
matrix elements, infrared and collinear cutoffs, angular ordering instead of gen- 
uine quantum interference, etc. Therefore, even if the hadron-level predictions 
of the model agree very well with experiment, as they often do, we cannot be 
sure that the hadronization effects are correctly estimated. 

In the past two or three years an alternative method has developed, based 
more directly on perturbation theory. The approach is a logical extension of 
the identification of suitable (infrared-safe or factorizable) observables within 
perturbation theory. One now asks which amongst these observables are ex- 
pected to be more sensitive to long-distance physics, and which less sensitive. 
This can be established by seeing whether low-momentum regions contribute 
much to the perturbative prediction, even when the overall momentum scale 
Q is high. Generally one finds that the low-momentum contribution is of the 
form A/QP, i.e., power-suppressed, modulo logarithms. Simply by examining 
Feynman graphs, one can compute the powers p for different observables, which 
is already useful. More speculatively, one can try to estimate the coefficients 
A, and to relate them to the power corrections observed experimentally. In the 
following Sections, I shall review recent work on both of these aspects of the 
problem. 

The approach outlined above has become known as the 'renormalon' or 'dis- 
persive' method for estimating power corrections. Sect. explains the origin of 
this terminology and its connection with low- momentum contributions. Then 
Sect, ra shows how the power corrections arise. Applications of the method to 
hadronic structure functions and event shape variables are reviewed in Sects. 

*As one of the engineers of this approach, I feel free to criticise it. 



and H respectively. Event shape variables are particularly 'bad' observables, in 
that they typically have the largest power corrections, with p = 1. In Sect. ^, 
I discuss the possibilities for constructing better shape variables, with p > 1. 
Finally Sect. H contains some conclusions and further speculations. 



2 Renormalons 

The usual starting point for a discussion of renormalons Q is the observation 
that the QCD perturbation series is not expected to converge. There may indeed 
already be signs of this in some of the quantities that have been computed to 
order a^, for example the Gross-Llewellyn Smith sum rule for deep inelastic 
neutrino scattering: 



dx{Fl^ + F^ 
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A known source of divergence is the class of diagrams in which a chain of 
(quark) loops is inserted into a gluon line as illustrated in Fig. 0. They lead to 
a factorial divergence of the form 



F^J2 ^»"s 
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(2) 



where b = —Nf/Gn for Nf flavours of quarks in the loops and the number p, 
normally an integer, depends on the observable F being computed. 




Figure 1: Diagram containing a renormalon chain. 

A crucial step in estimating the form of divergence of the full perturbation 
series at high orders is that of so-called 'naive non-Abclianization' [0, p|, which 



consists of making the replacement Nf — > Nf — 33/2 in the above equation, so 
that b becomes the first coefficient of the QCD /3-function, b = (33 — 2iV/)/127r. 
This takes account of gluon loops and other contributions that give rise to 
the running coupling. Surprisingly, the resulting estimates of the perturbative 
coefficients c„ seem to lie close to the true values even when Nf — (pure gauge 
theory), for quantities that have been studied to high orders by means of lattice 
perturbation theory Q. 

Assuming the approximate validity of the naive non-Abclianization proce- 
dure, the divergent series of contributions can be conveniently summarised by 
the formula g 

F{Q\ x)=J^ a^,,{fi') HfiVQ^ x) . (3) 

Here Q^ is the hard process scale, assumed large, and x represents dimcnsionless 
ratios of other hard scales, e.g. the Bjorkcn variable. As discussed more fully in 
[gl , ^^ is the variable of integration in a dispersion relation but can be regarded in 
some respects as an effective gluon mass-squared. The 'effective strong coupling' 
a^it is essentially equivalent to as at large values of fj^ but different from the 
perturbative expansion of as at low scales. This discrepancy might or might not 
depend on the observable F -we shall discuss later the consequences of assuming 
that a^H is universal. The 'characteristic function' J-" is found by evaluating the 
observable F to first order for a finite gluon mass ii. On dimensional grounds, for 
F dimcnsionless, this must be dependent on fi only via the quantity e = fi^/Q"^. 
T is then minus the logarithmic derivative of T with respect to e, i.e., with 
respect to p? . 

In general, ^ is a well-behaved function of e, vanishing like a power at small 
e and an inverse power at large e (modulo logarithms). Thus we expect the 
integrand in Eq. (O), considered as a function of the complex variable z = Ine, 
to be damped exponentially in both directions along the integration contour, 
which is the real axis (Fig. ||) . 
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Figure 2: Integration contour and convergence region in the Ine plane (e = 
To represent Eq. (0) perturbatively, we expand the effective coupling in terms 



ofas(g'): 

n 

a^lji^) = a.«(e^g') - ^ a„(z)K(Q')]" ■ (4) 

1 

This series has some finite region of convergence in the z-plane, iimited by the 
positions of the infrared singularities of a^a. If, for example, we use the one- loop 
expression arff(/i^) = l/61n(/x^/A^), then the region is a circle which touches 
the Landau pole at z = — ln((5^/A^). In this case a singularity lies on the 
contour of integration and the integral is not defined. In a more realistic model 
for Qeff, however, the infrared singularities will lie elsewhere in the complex z 
plane, as depicted in Fig. g, and the integral will be unambiguous order by order 
in asiQ"^)- Nevertheless the perturbation series in as(Q^) will diverge, because 
the contour extends outside the region of convergence, both to the right and 
to the left. Since T is exponentially damped in z on both sides, while a„(z) 
is power-behaved, with the power increasing with n, the divergence will be of 
factorial form. 

The divergence due to the segment of the contour to the right of the conver- 
gence region is called an ultraviolet renormalon. We can see directly from the 
integral representation (H) that this divergence is just an artefact of expanding 
in powers of as((5^). The ultraviolet contribution to the integral (say, that from 
z > 0, i.e. /x^ > Q^) is perfectly well defined in terms of the usual perturbative 
form of the running coupling acft(/i^) — as(/i^). It can be evaluated numeri- 
cally, or equivalently the corresponding series can be summed using a standard 
technique such as Borel summation. The divergence does not originate from 
the ultraviolet region, being simply a reflection of the limit on the region of 
convergence set by singularities far away in the infrared. 

On the other hand the divergence of the series due to the region on the left, 
called an infrared renormalon^ cannot be eliminated without information on the 
non-perturbative behaviour of a^it in the infrared region. It should be clear 
from Fig. that this ambiguity has nothing to do with whether a Landau pole 
is present on the axis H. It will exist as long as a^a has any kind of infrared 
singularities, which it surely must. Therefore the so-called infrared renormalon 
ambiguity can only be resolved by providing a physical model of the infrared 
behaviour of a^a. Once this is done, one can integrate unambiguously along the 
left-hand portion of the contour, since no physical singularities can lie on the 
z-axis (they would correspond to tachyons). 

3 Power corrections 

Assuming the validity of Eq. (pi), we can make useful deductions about the 
behaviour of F{Q^,x) at low Q , without any specific assumptions about the 
infrared behaviour of the effective coupling. We shall argue that F is of the 



form 

FiQ' , x) = fJ:] (g2 , x) + F^:} (Q2 , x) (5) 

where the first term, representing the standard n-th order perturbative predic- 
tion, has a basically logarithmic Q^-dependence while the second, which includes 
non-perturbative effects, is power-behaved. Note that the second term also re- 
quires a superscript n since it depends on where we truncate the perturbation 

(n) 

To derive Eq. (0) we split a^H into two parts, roughly described as pertur- 
bative and non-perturbative: 

a^.{f^') ^ ai:l{^l') + Sa^^^,\^,') . (6) 

(n) 

The superscript on a^Jt indicates that this quantity is obtained by expanding 
to n-th order in as{Q^)'. 



series F^,,^ 



c^i'lit^') = MQ^) + h[\n%+K]al{Q 



+ ... + o(\n^^ aUQ') (7) 

where the constant K depends on the choice of renormalization scheme. Once 
again it follows that the remainder Sa"; also depends on n. 

In Fig. |3| the solid curve shows a possible model for the behaviour of a^ff, 
in which it is assumed to vanish at very small scales, as suggested in ||7|]. Given 
the present level of understanding, it could equally well approach a non-zero 
value. The quantity Sal^/ is given by the difference between the solid curve and 

(n) 

the dot-dashed one representing a^Jt, which rises as n increases. Therefore we 
expect in general that 5ai2 will decrease with increasing n. 
Substituting in Eq. (g), we obtain Eq. (o) with 



i^iri = r K -^"-^ (m') ^(MVQ^ -) . (8) 
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We have assumed here that the approximation (ffl) to a^ff is sufficiently accurate 
at high fj? for Sal^f to be negligible above some 'infrared matching' scale /Xi 
( A <C A*i <C Q) . It follows then that the Q^ dependence of F^^l is determined 
by the behaviour of ^(e, x) for small values of e. For the general behaviour 



we have 



J^(e,a;)-/(x)eP/2in9e as e ^ (9) 

FU-^l ^'5«e«(M)M^ln^(^^j (10) 
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Figure 3: Possible behaviour of the effective coupling. Solid: «„«. Dashed: 
1-loop running as. Dot-dashed: ctplL for n = 2,3,4. 



Thus, as stated earlier, fJ"^ is power-behaved (modulo logs), with x-dependence 
given by that of .F(e, x) at small e and magnitude determined by a (log-)moment 
of the 'effective coupling modification' (5ai„ . 

All this depends only on the validity of the dispersive formula (ra) , irrespec- 
tive of whether 5al^{ varies from one observable F to another. Therefore one 
would expect the value of the power p to be given reliably by Eq. (^). Under 
the additional assumption that a^a, and hence (5aiJJ , is universal, one obtains 
testable relationships between power corrections to different observables with 
the same value of p. 

Over the last few years, this 'renormalon-inspired' or 'dispersive' approach 
to the estimation of power corrections has been applied to a wide variety of 
QCD observables Q,(|]-|l^. Many of the results have been presented in other 
talks at this meeting. Therefore I shall concentrate in the remainder of this talk 
on topics that illustrate some of the general remarks made above. 



4 Structure functions 

The need for 1/Q^ power corrections to hadronic structure functions, normally 
called higher- twist contributions, is well established and their dependence on 
the Bjorken x variable has been extracted from deep inelastic scattering (DIS) 
data. Typical results are shown in Figs. ^ and ^ |20: |l|. In Fig. ^ the higher- 
twist contribution is parametrized as D2{x)/Q^ time the leading- twist, whereas 
in Fig. ra it is represented simply as an additive term h{x)/Q'^ . 

Such contributions are predicted by the renormalon approach, since the char- 
acteristic functions for the evolution of (non-singlet) structure functions have 
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Figure 4: Higher-twist contribution to F2 (data, solid curve) and F3 (dashed 
curve) . 



the general form 



j=-{e,x) r^P{x)+f{x)e\ne + 



(11) 



where P{x) is the usual splitting function and the ellipsis represents terms 
that are less important as e -^ 0. The convolution of the (structure- function- 
dependent) coefficient function f{x) with the quark distribution then gives the 
predicted form of the leading 1/Q^ correction. 

The solid curve in Fig. ^ shows the prediction for the structure fmiction F2, 
with the relevant {p — 2, q = 1) moment of SaJ^t (for n = 2) in Eq. (Il^) as the 
only free parameter ||, |l^ |l^ . The good agreement with the data is due mainly 
to the characteristic 1/(1 — x) behaviour at large x, which is also predicted by 
a number of models of higher-twist contributions. 

A more stringent test is the prediction of the power correction to xF^ Q , 
which depends on the same moment of (Jaltf ■ The dashed curve shows the 
result, for the value of the moment obtained by fitting F2. It is interesting that 
the prediction is negative for all but the largest x values, in agreement with the 
need for a negative power correction to the GLS sum rule (0) . This is confirmed 
by a recent analysis |2l] of the data of the CCFR collaboration |g2[ : the power 
corrections extracted in LO, NLO and NNLO (i.e. n = 1,2,3) analyses are 
shown in Fig. ^, together with the same prediction for F3 as in Fig. ^. 

Fig. H also illustrates an important point discussed in Sect. ^, namely that 
the normalization of the power correction depends on the order of the pertur- 
bative analysis. It decreases in magnitude with increasing order, as more of the 
divergent renormalon series is included in the perturbative prediction Fpcrt ■ We 
can see from the GLS prediction, Eq. (|l|), that the perturbative coefficients are 
indeed increasingly negative, consistent with the presence of an infrared renor- 
malon with a negative coefficient, which tends to cancel the power correction. 




Figure 5: Higher-twist contribution to F^, extracted from various orders of 
perturbative analysis. Figure adapted from Ref. [21]. 



5 Event shape variables 

Power corrections to hadronic event shapes are generally larger than those in 
most other observables. This, together with the lack of complete next-to-next- 
to-leading order perturbative calculations of these quantities, has made it dif- 
ficult to use them for precise tests of QCD, in spite of the very high statistics 
available from LEPl. 

Usually the non-perturbative contributions to event shapes are estimated 
using Monte Carlo event generators [g3|, |2J] , or rather the hadronization models 
that are built into them. The difference between the parton- and hadron-level 
predictions of the generator is taken as the hadronization correction, with an 
error given by the discrepancy between the models. As discussed in Sect, m, 
the trouble with this is that the parton-level generators do not coincide with 
fixed-order perturbation theory. They represent an approximate summation of 



various contributions (including the renormalon graphs) up to rather high order. 
In addition, they contain cutoff parameters that can emulate non-perturbative 
effects. Thus the power corrections might be under- or over-estimated by this 
procedure, even though the generator might agree perfectly with the data at 
the hadron level. 

Some less model-dependent insight into power corrections to event shapes 
can be gained from the renormalon-inspired type of treatment discussed above. 
Consider as a prototype the thrust variable in e+e^ annihilation, 

r-max^|p,.n|/^|pj, (12) 

i i 

where the sum is over final-state hadron momenta p^ and the maximum is with 
respect to the direction of the unit vector n (the thrust axis). In lowest order the 
difference 1 — T arises from a qqg final state. Making a Sudakov decomposition 
of the gluon momentum fc, 

k = ap + I3p+ k±^ (13) 

where p^ = P^ — P ■ k±^ = p ■ ki_ =0, k\ = —k\ and 2p ■ p = Q^, the cm. 
energy-squared, we have 

t = l-T = min{a, (3} (14) 

in the soft region. Hence the characteristic function for the mean value of this 
quantity is 

TT J a (3 
■ min{a, (3} S{af3Q'^ - k\ - ^i^) 

(15) 



2Cf f dk\ 



t:Q 



ki+fi^ 



^-^^. (16) 

TT Q 



from which we find 



It follows from Eqs. (g) and (^^ that the quantity (t) is expected to receive a 
1/Q power correction, in contrast to the l/Q^ corrections to other quantities 
like structure and fragmentation |15| functions. A similar result is obtained for 
the mean values of several other e+e^ shape variables |g, |l^, |^, and for the 
corresponding quantities in deep inelastic final states [|l6| . 

There is by now good experimental evidence for 1/Q corrections to the mean 
values of several event shape variables in both e+e^ annihilation |2^, [2^ and DIS 
[ETI . Furthermore the magnitudes of the corrections are in most cases consistent 
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(at about the 20% level) with a universal value of the relevant moment of 6al2 
in Eq. ([lO|), namely 

^^^Sa^T'\^^')^lGeY. (17) 



M 



Turning from the mean values to the distributions of event shape variables, 
a simple result is obtained when the dominant low-momentum contributions 
exponentiate, that is, when they correspond to essentially uncorrelated multiple 
soft gluon emission. This is the case for the thrust distribution, which may be 
expressed as a Laplace transform pqj 

F{t) = --^ = / di^ exp[i^t + R{v, Q)] (18) 

a at Jq 

where the contour C runs parallel to the imaginary axis, to the right of all 
singularities of the integrand, and f = 1 — T as usual. The 'radiator' R has the 
following form in lowest order: 



>^,Q) = 


2Cf fd^I^ . 2, 

TT J M 




. / ^ (e— - 



(19) 

To compute the effect of a non-perturbative contribution to a^a, confined to the 
region fj?' < fi^ <^ Q^, we note that, since v is conjugate to i, the exponential 
can safely be expanded to first order as long as t :$> fi^/Q. The resulting change 
in the radiator is 



2Cf f dfl'^ r (n)/ 2^ Ai ^ 



5R^--^v / '^5a^2\^,^)^ ^ -v^ (20) 



which corresponds simply to a shift in the t-distribution ^, |l^ , by an amount 
AjQ equal to the non-perturbative shift in (i) computed above: 

F{t)^F^^,,{t-AIQ) (21) 

where -Fpoit is the perturbative prediction. This agrees well with the e+e^ data 

There is some ambiguity in the calculations outlined above: we assumed an 
effective gluon mass [i but used the massless forms of the shape variable and 
lowest-order matrix element, when expressed in terms of the Sudakov variables a 
and /3, as was done in B. In reality the 'massive' gluon fragments into massless 
partons, which eventually form hadrons, and the fragmentation products may 
have a thrust value different from that of their parent |13] . This will not change 
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the l/Q form of the power correction but it would be expected to modify its 
coefficient. 

These issues have been clarified in a recent paper from the Milan group p9 , 
in which the coefficient of the 1/Q correction is evaluated to two- loop order in 
the soft limit. A substantial increase, relative to the above naive treatment, is 
obtained - a factor of 1.8 for three active flavours. Interestingly, however, this 
enhancement factor itself appears to be universal across a wide range of shape 
variables, in both e+e^ and DIS final states ^^- Therefore the experimental 
analyses |g5|, |2^, |2^ of event shapes assuming universality are not invalidated, 
but the fitted values of the moments of (5aL" , such as Eq. ([l7|), need to be 
rescaled by the inverse of the Milan factor. 

6 Better shape variables 

As we have seen, the power corrections to event shapes are normally of order 
l/Q. For determinations of ag, it would be an advantage to know if there are 
any shape variables that are free of such corrections, preferably up to order 
l/Q\ 

One set of variables that has been suggested is the higher moments of shape 
variables. A straightforward approach along the lines of Sect. S suggests that if 
the mean value has a correction of order 1/Q then the p-th moment will only 
receive a non-perturbative contribution of order 1/Q^. 

Unfortunately the situation is not so simple, since we would also like to 
safeguard against corrections of order as{Q^)/Q ^ 1/Q^- Consider again for 
example the variable t = 1 — T. We know that the power correction to the NLO 
perturbative prediction for (t) is A/Q where A ~ 1 GeV. Furthermore we saw 
that the main non-perturbative effect on the thrust distribution F{t) is a simple 
shift, Eq. (pi|). It follows that for the mean value of fP we have 

(tP) - j dttPF{t)^ [ dt{t + A/QYF^Ut) 



P^t^p-U : ,nf ^ 



= (^')p- + ^(^'"')p-+0(^Qjj (22) 

In practice, the perturbative predictions for {t^) decrease rapidly with p, and 
therefore the second term is relatively important for p > 1, in spite of its sup- 
pression by as{Q^)- One finds 

(t)pert = 0.335as + 1.03a^ -I- . . . , 
(i^)port = O.OSOas + 0.15a^ + . . . . (23) 

and so, at LEP energies and below, the power corrections to {t) and (P) are 
actually comparable, relative to the perturbative predictions. In addition, the 
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perturbation series for (i^) and higher moments show increasingly poor conver- 
gence. 

A better idea is to consider moments about the mean value, such as the 
thrust variance a^ = (t^) — (i)^. Then we see from Eq. ( |2^ ) that the as{Q^)/Q 
corrections should cancel out. Furthermore the perturbative prediction, 

o"?, port = O.OaOas + OMaj + ... , (24) 

shows better convergence than that for (i^) - better than (t), in fact. Nev- 
ertheless, this is a quantity an order of magnitude smaller than (t), and so 
its measurement will require that experimental sources of error are well under 
control. 

An interesting open question is whether the 1/Q and as{Q'^)/Q corrections 
to a^ remain cancelled after taking into account the Milan enhancement factor, 
to be extracted from a two-loop analysis. 

Finally we should consider a shape variable which is not expected to have 
1/Q or as{Q^)/Q corrections to its first moment, namely the three-jet resolution 
2/3. This is the value of the Durham jet measure y<.ut at which three jets are just 
resolved in the final state. Recalling that the Durham criterion for resolving 
two partons i and j is yij > y^^t , where [^ 

y,j = 2inm{El E^}{1 - cos0,j)/Q^ , (25) 

we find that in the soft limit we obtain in place of Eq. ( [l4| ) 

y3 = {a + f3)mm{a,(3}, (26) 

giving instead of Eq. (fi^ ) 

for (ys). Thus, modulo possible logarithms, the expected leading correction to 
this quantity is 0{1/Q^). Once again, it would be interesting to know whether 
the Milan enhancement can somehow generate a 1/Q correction. Experimen- 
tally, no significant power correction to (2/3) has yet been found |2(| , which is 
certainly consistent with the absence of a 1/Q term. 

7 Conclusions 

The renormalon-inspired or dispersive method for estimating power corrections 
to perturbative calculations of QCD observables has proved to be a useful phe- 
nomenological tool. At the very least, it warns us of those quantities and phase- 
space regions in which the perturbative calculation is dangerously sensitive to 
low-momentum contributions. If we wish, we can simply avoid those regions in 
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comparing theory with experiment. From this viewpoint, the renormalon ambi- 
guity is analogous to the renormahzation-scale dependence of a next-to-leading 
order perturbative prediction: if it is large, we would be unwise to trust the pre- 
diction as a reliable test of QCD or as a good way of measuring as ■ The ideal 
observable for these purposes has both a small scale dependence and a small 
renormalon ambiguity - in particular, a large value of the power p. These are 
somewhat independent requirements, and therefore a renormalon analysis pro- 
vides useful additional information, even if we do not accept it as a quantitative 
estimate of power corrections. 

As far as I am aware, the powers p for the corrections observed experimen- 
tally agree with those predicted in all cases studied so far. This provides the 
motivation for going further and comparing the measured and computed co- 
efficients of power corrections, assuming a universal low-energy form for the 
effective strong coupling a^ti- The initial results look encouraging to me, al- 
though many issues remain unclear. As emphasised in Sect. H, the coefficient 
must depend on the order of the perturbative prediction that is being sub- 
tracted. Presumably the optimal comparison would be with the perturbation 
series truncated at its smallest term, but we remain far from that situation. 
In the case of event shape variables, the renormalization of the coefficient by 
two-loop effects has been shown to be large but appears to remain universal 
for a range of observables. A general argument can be made that higher loop 
corrections should be small ||2^, but explicit calculations would be reassuring. 

Another relevant question is whether we should expect the strong coupling 
itself to exhibit power corrections, possibly as large as 1/Q^ p2[ . There have 
been recent suggestions of such terms in lattice calculations [p3[. A systematic 
perturbative treatment of the effective strong coupling beyond one-loop (i.e. 
one-renormalon-chain) order seems essential for this purpose, but is still lacking. 
Approaches such as that of gj] may be fruitful. Such a treatment would also be 
valuable for the study of power corrections to gluon-dominated quantities such 
as hadronic jet shapes |35[| and small- a; structure functions |36|1 . 
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